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We consider the kick (tilt)-induced mobility of two-dimensional (2D) fundamental dissipative 
solitons in models of bulk lasing media based on the 2D complex Ginzburg-Landau (CGL) equation 
including a spatially periodic potential (transverse grating). The depinning threshold, which depends 
on the orientation of the kick, is identified by means of systematic simulations, and estimated by 
means of an analytical approximation. Various pattern-formation scenarios are found above the 
threshold. Most typically, the soliton, hopping between potential cells, leaves arrayed patterns of 
different sizes in its wake. In the single-pass-amplifier setup, this effect may be used as a mechanism 
for the selective pattern formation controlled by the tilt of the input beam. Freely moving solitons 
feature two distinct values of the established velocity. Elastic and inelastic collisions between free 
solitons and pinned arrayed patterns are studied too. 

PACS numbers: 42.65.Tg,42.65.Sf,47.20.Ky 
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I. INTRODUCTION 

A well-known fact is that the formation of stable dissipative solitons — most typically, in lasing media 0, 0] and 
plasmonic cavities Q — relies upon the simultaneous balance of competing conservative and dissipative effects in the 
system, i.e., respectively, the diffraction and self- focusing nonlinearity, and linear and nonlinear loss and gain Q. The 
generic model describing media where stable dissipative solitons emerge via this mechanism is based on the complex 
Ginzburg-Landau (CGL) equations with the cubic-quintic (CQ) combination of gain and loss terms, which act on top 
of the linear loss Q. In addition to modeling the laser-physics and plasmonic settings, the CGL equations, including 
their CQ variety, serve as relevant models in many other areas, well-known examples being Bose-Einstein condensates 
in open systems (such as condensates of quasi- particles in solid-state media) reaction-diffusion systems @, and 
superconductivity Q. Thus, the CGL equations constitute a class of universal models for the description of nonlinear 
waves and pattern formation in dissipative media @. 

The CGL equation with the CQ nonlinearity was originally postulated by Petviashvili and Sergeev Q as a model 
admitting stable localized two-dimensional (2D) patterns. Subsequently, systems of this type were derived or intro- 
duced phenomenologically in many physical settings, and a great deal of ID and 2D localized solutions, i.e., dissipative 
solitons, have been studied in detail in them [Ioj-jl3|. 

A 2D model of laser cavities with an internal transverse grating, based on the CQ-CGL equation supplemented 
by a spatially periodic (lattice) potential, which represents the grating, was introduced in Ref. [ljj]. Note that the 
currently available laser- writing technology makes it possible to fabricate permanent gratings in bulk media (l5j . In 
addition, in photorefractive crystals virtual photonic lattices may be induced by pairs of pump laser beams with 
the ordinary polarization, which illuminate the sample in the directions of x and y, while the probe beam with the 
extraordinary polarization is launched along axis z |16j . In fact, the laser cavity equipped with the grating may be 
considered as a photonic crystal built in the active medium. Periodic potentials are also known in passive optical 
systems, driven by external laser beams and operating in the temporal domain, unlike the spatial-domain dynamics 
of the active systems. In such systems, effective lattices may be induced by spatial modulation of the pump beam 
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A notable fact reported in Ref. [14| is that localized vortices, built as sets of four peaks pinned to the periodic 
potential, may be stable without the presence of the diffusion term in the CGL laser model, which is necessary for 
the stabilization of dissipative vortex solitons in uniform media, see, e.g., Ref. [Tl| . but is unphysical for waveguiding 
models (the diffusion term is relevant in models describing light trapped in a cavity, where the evolutional variable is 
time, rather than the propagation distance [IH). In subsequent works, stable fundamental and vortical solitons in 2D 
f20j and 3D [2l| CGL models with trapping potentials were studied in detail. Spatiotcmporal dissipative solitons in 
the CQ-CGL model of 3D laser cavities including the transverse grating were investigated too [2l[. Both fundamental 
and vortical solitons were found in a numerical form as attractors in the latter model, and their stability against 
strong random perturbations was tested by direct simulations. 

While the stability of various 2D localized patterns has been studied thoroughly in the framework of the CQ-CGL 
equations with the transverse lattice potential used as the stabilizing factor, a challenging problem is mobility of such 
2D dissipative solitons under the action of a kick applied across the underlying lattice. 

Actually, the action of the kick in this context implies the application of a tilt to the beam. It should be noted that 
the CGL equation models single-pass optical amplifiers, as well as laser cavities. In the latter case the evolution is 
considered in the temporal domain, while in the former the evolution variable is the propagation distance. In order 
to build the required initial data in a laser cavity, the tilt should be applied very quickly. This may be achieved 
by means of an optically induced grating or mirrors, using some fast device. On the other hand, in the single-pass 
amplifier the input is an incident beam, the tilt being a mere misalignment between the beam and the amplifier's axis. 
Below, we present the analysis in terms of the amplifier's setup, which is more straightforward for the experimental 
implementation. 

Thus, we assume that an external source produces a beam, which is shaped into the fundamental-soliton mode of the 
amplifier by means of an adequate setup. Then, the direction of this input beam is slightly tilted with respect to the 
amplifier's axis, allowing the transverse part of the beam's momentum to acts as the kick applied to the fundamental 
spatial soliton. 

Furthermore, we demonstrate that the effective hopping motion of the kicked soliton through cells of the periodic 
potential can be used for controlled creation of various patterns filling these cells (or a part of them). 

Thus, the main objective of this work is to study the mobility of the 2D fundamental solitons, and scenarios of the 
pattern formation by kicked ones, in the framework of the CQ-CGL models with the lattice potential. The model 
is formulated in Section II, which also presents an analytical approximation, that, using the concept of the Peierls- 
Nabarro (PN) barrier, makes it possible to predict, with a reasonable accuracy, the minimum (threshold) strength of 
the kick necessary for depinning of the quiescent soliton trapped by the lattice. The main numerical results for the 
mobility of the kicked soliton, and various scenarios of the pattern creation in the wake of the soliton hopping between 
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cells of the potential lattice, are reported in Section III, while Section IV deals with collisions between a freely moving 
soliton and a standing structure created and left by it, in the case of periodic boundary conditions (which correspond 
to a pipe-shaped amplifier, i.e., one in the form of a hollow cylinder). In particular, elastic collisions provide for an 
example of a soliton Newton's cradle. The paper is concluded by Section V. 



II. THE MODEL AND ANALYTICAL APPROXIMATIONS 



A. The Ginzburg-Landau equation 

Following Refs. [l4| and [2(|, the scaled CQ-CGL equation for the evolution of the amplitude of the electromagnetic 
field, u(X,Y, Z), in two dimensions with transverse coordinates R = (X, Y), along the propagation direction, Z, is 
written as 



da 
dZ 



-S + tV{X,Y) + -Vl + {i + e)\u\ 2 - {w + ^)\u\ 4 



(1) 



where the paraxial diffraction is represented by Vj_ = d 2 /dX 2 + d 2 jdY 2 , real coefficients S, e, and /i account for 
the linear loss, cubic gain, and quintic loss, respectively, and coefficient v > accounts for the saturation of the Kerr 
nonlincarity. The transverse grating is represented by the periodic potential, 

V(X, Y) = -Vo [cos(2X) + cos(2Y)] , (2) 

of depth 2Vo, with the period scaled to be tt. Localized modes produced by Eq. ([T]), which physically correspond to 
light beams self-trapped in the (X, Y) plane, are characterized by the total power, 

P = J J \u(X,Y,Z)\ 2 dXdY. (3) 

In the simulations, Eq. ([1]) was solved by means of the standard fourth-order Runge-Kutta scheme in the in- 
direction, and five-point finite-difference approximation for the transverse Laplacian. As specified below, we used 
periodic boundary conditions. The integration domain correspond to an (X, Y) matrix of 256 x 256 grid points 
covering the area of \X\ , \Y\ < 22. Generic results for the mobility of fundamental dissipative solitons can be 
adequately represented by fixing the following set of parameters: 

S = 0.4, e = 1.85, fj, = l,v = 0.1, V = 1, (4) 

for which the quiescent fundamental soliton is stable. 



B. The description of tilted beams 



In simulations of Eq. ([T]), the kick (i.e., tilt) was applied to the self-trapped beam, as usual, by multiplying the 
respective steady state, uq, by exp (ikoR), with the vectorial strength of the kick defined as 

ko = (ko cos 8, ko sin 8) , (5) 

where the square-lattice symmetry of potential ([2]) makes it sufficient to confine the orientation angle, 8, to interval 
< 8 < 7r/4. In terms of the amplifier setup, a small angle ip between the carrier wave vector Ko = (K x , K y , K z ) of 
the beam (in physical units) and the Z-axis, gives rise to the transverse component ko which corresponds to the kick. 
The paraxial approximation implies that 

if » sjKl+iq/K z « 1, (6) 

Before studying effects induced by the kick, it is relevant to explain the corresponding physical setting in more detail, 
making it sure that the tilted beams remain within the confines of the paraxial description. 

Equation (JXJ) can be derived from the underlying wave equation by means of the standard slowly varying envelope 
approximation (SVEA) [l[. For this purpose, electric field E (in its scalar form) is split into a slowly varying amplitude 
and the rapidly oscillating carrier either as 



E = A{x, y,z- vt y( K ^+ K vy+ K ^-^) + ex., 



(7) 
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where (x,y,z,t) are the coordinates and time in physical units, v is the group velocity, and c.c. stands the complex 
conjugate, or, alternatively, as 

E = A(x, Vl z~ vty( K * z - ut) + c.c. (8) 

The difference is that that carrier wave is oblique in Eq. ([7]) , while in Eq. (jSJ it is always defined as the straight one. 
Obviously, A(x,y, z — vt) = A(x, y, z — vt)e l ^ KxX+KvV ^ , and the two forms are fully equivalent within the framework of 
the SVEA if the oscillations due to the term e % (. K * x + K vv) are not essentially faster than the transverse variations of the 
beam described by amplitude A(x, y,z — vt). Thus, the SVEA can be fixed in the form of Eq. ©, with K z = 27m/A, 
where n the linear index of the medium. 

Further, to proceed to normalized equation ([!}, we set {X, Y} = w^ 1 {x, y}, where w is a characteristic scale of the 
beam's width — say, w ~ 10 (im for narrow beams, if the underlying wavelength (in vacuum) is A ~ l ^m. Accordingly, 
the period of grating ^ is Ttw in the physical units. Note that w ~ 10 fim corresponds to the period nw ~ 30 firm, 
which is relevant estimate (gratings with the period on this order of magnitude can be readily manufactured). Then, 
the scaled propagation distance is Z = z/(K z w) = (\/2irnw) z, and the scale wave amplitude is 

. = ^Re( x (3))^ 

where x*- 3 ' is the third-order nonlinear susceptibility. Further, the respective rescaling of the wave vector components 
is given by {k x ,k y } = w {K x , K y }. Thus, the deviation angle can be estimated as tamp = \J K 2 + Ky I = 

(X/2Trnw) ^Jk 2 + k 2 i and, for generic tilted modes, with yjk"* + ~ 1 in the scaled notation, condition ^ for the 

validity of the paraxial approximation amounts to A 2imw, which is nothing but the standard paraxial assumption. 

With the above-mentioned typical values, A ~ 1 /mi and w ~ 10 fim, along with n « 1.5, the above estimate yields 
ip ~ 0.01 (in radians). Then, the minimum propagation distance relevant to the experiment, z ~ 1 cm, corresponds to 
the transverse deviation of the tilted beam Ax ~ 100 /im, which can be easily detected and employed in applications. 



C. Analytical estimates 



The first characteristic of kink-induccd effects is the threshold value, (fco) thr , such that the soliton remains pinned at 
fco < (fco)thr an d escapes at fco > (fco) thr .To develop an analytical approximation which aims to predict the threshold, 
one can, at the lowest order, drop the loss and gain terms, as well as the lattice potential, in Eq. (TT]). The corresponding 
2D nonlinear Schrodinger equation gives rise to the commonly known family of Townes soliton, which share a single 
value of the total power, Pt ~ 5.85 [22[. The family may be approximated by the isotropic Gaussian ansatz with 
arbitrary amplitude A, 

u{z,R= VX 2 + Y 2 ^) = A cxp (ibZ -T^rR 2 ^), (9) 

and propagation constant b = A 2 /4 [23[. Then, taking into regard the loss and gain terms as perturbations, one can 
predict the equilibrium value of the amplitude from the power-balance equation, 

8-P + lil [ \u(X,Y)\ 6 dXdY = e [ I \u{X,Y)\ A dXdY. (10) 



(11) 



The substitution of approximation into Eq. (|10p leads to a quadratic equation for A 2 , with roots 

2 _ 3e ± y/3 (3e 2 - 16 fiS) 



the larger one corresponding to a stable dissipative soliton (cf. a similar analysis for the CQ-CGL model in ID [24|). 

Proceeding to the kicked soliton, the threshold magnitude of the kick for the depinning, (fc ) tllr , can be estimated 
from the comparison of the PN potential barrier, t/pN, and the kinetic energy of the kicked soliton, i^kin- The Galilean 
invariance of Eq. (JTJ) [25[ (in the absence of the lattice potential) implies that kick ([5]) gives rise to "velocity" ko, so 
that the solution will become a function of R — ko-Z = R T, instead of R, the corresponding kinetic energy being 

E kin = (1/2)P*§ (12) 
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(P plays the role of the effective mass of the soliton). Further, assuming that the kicked soliton moves in the direction 
of 9 [see Eq. (|3J|], the effective energy of the interaction of the soliton, taken as per approximation ©, with lattice 
potential ((2|), treated as another perturbation, is 

E pot {T) = -V [ J [cos(2X) + cos(2y)] \u(X - Tcos 9,Y - Tsin 9) | 2 dXdY 

= -VoPcxp ( -— £r | [cos (2T cos 9) + cos (2T sin 0)1 , (13) 

where T is the shift of the soliton from X = Y = in the direction of 9. As follows from this expression, the PN 
barrier, i.e., the difference between the largest and smallest values of the potential energy, is estimated as 

E PN = V Pexp(--^j A(0), (14) 

where A(0) is the difference between the maximum and minimum of function cos (2T cos 6) + cos (2T sin (9). Obviously, 
A(0) = 2 and A (7r/4) = 4. For intermediate values of 9, it may be approximated by the difference of the values of 
the function between points T = and T = ir/ (2cos0), i.e., 

A(0) w 3 - cos (tt tan 9) . (15) 

Finally, the threshold value of the kick is determined by the depinning condition, -Ekm = -Epn, i.e., 

(*o)thr = A/2y o [3-cos(7rtan0)]exp (-^2) ■ ( 16 ) 
This prediction is compared with numerical results below. 



III. NUMERICAL RESULTS: THE MOBILITY AND PATTERN FORMATION 



The stable fundamental soliton constructed in the model based on Eqs. (JTJ and ([2]) at parameter values (j4} is 
shown in Fig. [TJ For these parameters, analytical prediction (fTTj) yields the amplitude of the stable soliton A i=s 1.496 , 
which is quite close to the amplitude of the numerically found solution in Fig. [TJ A *=s 1.479 , which implies that the 
isotropic Gaussian (J9j> is quite appropriate as the ansatz for the description of static properties of the solitons. 



>H0 




(a)Amplitudc (b)Profilc 

FIG. 1. (Color online) The stable fundamental soliton: (a) the contour plot of the local amplitude, \u(X,Y) |; (b) the cross- 
section profile of |u(-X^| at Y = 0. The soliton emphasizes the potential : V(X, Y) = — [cos (2X) + sin (2X)] 



A. The formation of arrayed soliton patterns 

First we consider the solitons kicked with = 0, i.e., along bonds of the lattice, see Eq. ([5]). Below the threshold 
value of the kick's strength, whose numerically found value is 

(fco)thr (0 = 0)- 1-6865, (17) 

the kicked soliton exhibits damped oscillations, remaining trapped near a local minimum of the lattice potential, as 
shown in Fig. [21 Originally (at < z < 8 in Fig. [5]), the total power ^ increases, and then it drops to the initial 
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value, P rj 3.2. As a result of the kick, a portion of the wave field passes the potential barrier and penetrates into the 
adjacent lattice cell, but, at k < (fco) thr , the power carried by the penetrating field is not sufficient to create a new 
soliton, and is eventually absorbed by the medium. 

On the other hand, analytical prediction (fTB"| yields (fco) thr (0 = 0) rj f .32. A relative discrepancy ~ 20% with 
numerical value (|17[) is explained by the fact that, near the depinning threshold, the moving soliton suffers appreciable 
deformation, while the analytical approach assumed the fixed shape (J9j) , and it did not take into account energy losses 
(the latter factor makes the actual threshold somewhat higher). In other cases considered below, see Eqs. (fT5)l and 
(|I9[) . the analytical predictions for (&o) thr is also ~ 20% smaller than their numerically found counterparts 




z z 



(a) (b) 



FIG. 2. (Color online) (a): The evolution of the local amplitude, \u(X, Z)\, in the cross section, Y = 0, of the fundamental 
soliton kicked with the under-threshold strength, ko = 1.61 at 6 — 0. (b) The evolution of the total power in this case. The 
horizontal line designates the power of the quiescent fundamental soliton. 

If the kick is sufficiently strong, k > (fco) thr , the portion of the wave field passing the potential barrier has enough 
power to create a new dissipativc soliton in the adjacent cell. The emerging secondary soliton may either stay in its 
cell, or keep moving through the lattice. 





(a)Z = 2.65 



(b)Z = 4.70 




(c)Z = 8.51 
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(d)Z = 14.00 



FIG. 3. (Color online) The evolution of the local field amplitude, \u(X, Y)\, corresponding to the kicked soliton, for ko = 1.6878 
and 8 — 0. The field distributions are displayed at different values of propagation distance Z. 

Figure [3] demonstrates the creation of two new solitons at ko = 1.6878, which slightly exceeds the threshold value 
(|17jl . This figure represents a generic dynamical scenario, that can be summarized as follows: 

• the initial soliton (or a part of it) passes the potential barrier and gets into the adjacent (second) cell; 
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• it then stays for some time in that cell; 

• if the initial kick is not strong enough, the secondary soliton permanently stays at this location; 

• if the kick is harder, the soliton again passes the potential barrier, getting into the third cell, and may continue 
to move through the grating; 

• a portion of the wave field of the secondary soliton stays in the second cell and grows into a full soliton in this 
cell; 

• if the kick is not sufficient to continue the filling of farther cells, oscillations of all the persistent solitons relax. 

The present case is further illustrated in Fig. 0]by the plot for the evolution of the total power, which shows that 
P attains the first maximum at Z = 15.89, and then oscillates. Every minimum correspond to the collision between 
the two solitons. For the sake of comparison, four horizontal lines in the figure mark the powers corresponding to the 
single stable soliton (P so \ w 3.15), multiplied, respectively, by 1, 2, 3, or 4. 



4 
3 
2 
1 



"0 50 100 150 200 250 300 
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FIG. 4. (Color online) The evolution of the total power of the pattern in the case displayed in Fig. [3] with fco = 1.6878 and 
9 — 0. The horizontal lines show the powers for the sets of 1, 2, 3 and 4 quiescent fundamental solitons. 

It has been found that solitons can duplicate several times, thus forming extended patterns in the form of soliton 
arrays. The increase of the kick's magnitude leads to the decrease of the number of the solitons forming this pattern, 
as the soliton moves faster and does not spend enough time in each cell to create a new soliton trapped in it. In 
particular, Fig. [5] demonstrates that only one additional soliton is generated at fco = 1.6872, both solitons remaining 
pinned (note that this value is smaller than fco = 1.6878 appertaining to Fig. [3]). Further, at fco = 2.082 the soliton 
performs unhindered motion, without leaving any stable pattern in its wake (not shown here in detail). 

On the other hand, the smaller kick can initiate the creation of an arrayed pattern. This outcome of the evolution 
is shown in Fig. [6j where the array of five solitons is created, starting with the soliton initially kicked by fco = 1.694, 
in addition to which a free soliton keeps moving as a quasi- particle (cf. Ref . [263 ) > until it collides with the array from 
the opposite direction, due to the periodic boundary conditions along X, and is subsequently absorbed by the array 
(the collision is displayed in panel |6jf), where an additional soliton is observed at X < 0). This dynamical scenario 
is additionally illustrated below by Fig. [TTJ 

The emerging array remains in an excited state, featuring localized density waves running across it, as shown, on a 
much longer scale of Z, by means of the cross-section picture in Fig. [7])- It is worthy to note that the wave is reflected 
from the last pinned soliton. Such localized density perturbations propagating through a chain of pinned solitons are 
similar to the so-called " super-fluxons" , which were investigated experimentally and theoretically in arrays of fluxons 
(topological solitons, representing magnetic-flux quanta) pinned in a long Josephson junction with a periodic lattice 
of local defects [27j , as well as in an array of mutually repelling solitons forming a Newton's cradle in a two-component 
model of binary Bose- Einstein condensates [28| . 




B. The dependence of the outcome of the evolution on the strength of the initial kick 

Results of the systematic analysis of the model are summarized in Fig. where the number of solitons in the 
stable arrayed patterns established by the end of the simulation, is plotted versus the initial kick fco for 6 = 0, where 
intervals of the values of fco corresponding to constant numbers of the solitons are adduced. 

In case the free soliton collides with the quiescent array after performing the round trip in the system with the 
periodic boundary conditions, see the example above in Fig. [51 and an additional one (for four solitons) in Figs. O 
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FIG. 5. (Color online) The same as in Fig. S but for k = 1.6872 and = 0. 



and 1101 the number of solitons was counted just before the first such collision. Otherwise, the number was recorded 
after any motion in the system would cease. 

The results demonstrate that the number of solitons in the established patterns rapidly increases from 1 at ko < 
(fco)thr ~ 1-6865 [see Eq. (|T7|) ] to a maximum of five solitons, plus a sixth freely moving one, at ko = 1.6927. New 
solitons add to the established pattern according to the scenario outlined by means of the bullet items in the previous 
subsection. At ko > 1.6927, the soliton number decreases by steps with an increasing length of the corresponding 
intervals of the kick's strength, see Fig. [8] 

As mentioned above, the largest number of six solitons is attained at 1.6927 < ko < 1.6942. In addition to Figs. 



[6] and [3 this situation is illustrated by Fig. 11(a) I, where the averaged total power is P — 23, see panel fTTTb) . 



The horizontal reference lines in this figure show the power levels corresponding to a single quiescent soliton (recall 
Pso\ ~ 3.15), multiplied by factors from 1 to 7, which demonstrates that the total power of the six-soliton complex 
exceeds the seven-fold power of the single soliton. This is due to the fact that the energy of the moving soliton is 
roughly twice that of a soliton at rest. 



At ko > 2.082, the kicked soliton moves freely across the simulation domain, see Fig. 12(a) In this case, the figure 
demonstrates that the soliton's velocity increases, approaching a certain limit value. The computation of the velocity 
was performed by means of the Lagrange's interpolation of the numerical data to accurately identify the soliton's 
center. To display the results, small-scale oscillations of the velocity of the soliton passing the periodic potential 
(obviously, the velocity is largest and smallest when the solitons traverses the bottom and top points of the potential, 
respectively) have been smoothed down by averaging the dependence over about fifteen periods. 



C. The evolution initiated by an oblique kick 

The application of the kick under an angle to the lattice, i.e., with 6^0 [see Eq. was considered too. The 
results are summarized in Table U for 6 = 7r/4 (the kick oriented along the diagonal), and, below, in Table [TTI for 
9 = tt/8. 

For 6 = 7r/4, the initial soliton remains pinned at 

ko < (fco) thr (0 = tt/4) = 2.166, (18) 

and it is destroyed at 2.25 < ko < 4.337. At ko > 4.337, the single soliton survives, moving freely along the diagonal 
direction. Thus, the final number of the solitons in this case is 1 or (see Table [J), and no new solitons are generated. 
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FIG. 6. (Color online) The same as in Figs. [3] and [5] but for ko = 1.694. In this case, the kicked soliton eventually creates an 
arrayed set of five solitons (in panel (f), the additional freely moving (sixth) soliton hits the array from the opposite direction, 
completing its round trip in the system). 



Number of solitons 


Range of ko 


1 


k € [0, 2.355] 





k £ [2.36, 4.337] 


1 


k e [4.563, 10] 



TABLE I. The number of solitons in the established pattern versus the kick's strength, fco, at 8 = n/4. 



As concerns the comparison with the analytical prediction (fT6]l , it yields (fco)thr = "V^) ~ 1-87, which, as well as 
in the case of 9 = 0, is somewhat lower than its numerical counterpart. 

The kick with much larger values of fco (one or two orders of magnitude higher than in Table [TJ) causes the generation 
of dark-soliton structures, supported by a nonzero background filling the entire domain (the total power may then 
exceed that of the single soliton by a factor ~ 1000). Inspection of Fig. 13 (a) | reveals stable holes in the continuous 
background, whose centers coincide with phase singularities (see Fig. 
by the finite background. 

The creation of new solitons is possible in the case of the oblique kick with 6 = tt/8. In this case, the new solitons 



13(b) ). These holes represent vortices supported 
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FIG. 7. (Color online) Evolution of the local amplitude, \u(X,Z)\, in the cross-section plane, Y — 0, in the same case as in 
Fig. M(ko = 1-694 and 9 = 0). 
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FIG. 8. (Color online)The number of solitons in the established pattern versus the kick's strength, fco, at 6 = 0. In narrow 
gaps between intervals presented in the graphic, the number cannot be defined exactly, as it changes there by 1 



may be oriented along either axis, X or Y, as indicated in Table IH (the total count includes the obliquely moving 
originally kicked soliton). 



Number of solitons 


Range of ko 


Total 


Along X 


Along Y 


1 








fco e [0, 1.816] 


3 


2 





fco = 1-974 


2 


1 





k = 2.1 


1 








fc 6 [2.224, 4.569] 











fco e [4.816, 5.804] 


1 








fco e [6.05, oo) 



TABLE II. The same as in Tabled] but for the oblique kick oriented under angle 6 — n/8. 
Further, for the kick applied at angle 9 = 7r/8, the simulations demonstrate that the kicked soliton remains pinned 

at 

ko < (fc o ) thr (0 = 7r/8) = 1.816, (19) 

while analytical approximation (|16|) for the same case predicts (fco) t hr = 7r /^ > ) ~ 1-54. The creation of new solitons 
occurs above the threshold, similar to the case of 9 = and in contrast with 9 = 7r/4. The largest three-soliton 
pattern is created at fco = 1.974. It is composed of two oscillating pinned solitons and a freely moving one, as shown 
in Fig. 1141 At fco = 2.1, the dynamics again amounts to the motion of a single soliton. Further, the simulations 
demonstrate that, in all the cases of the free motion of the single soliton, it runs strictly along the X-axis, despite the 
fact that the initial kick was oblique. 
A harder kick, 



4.816 < fc < 5.804, 



(20) 
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FIG. 9. (Color online) The same as in Fig. [6l but for ho = 1.705. In this case, a free soliton splits away from the quiescent 
array, then hits it from the other side, and gets absorbed by it. 




50 100 150 200 250 300 



FIG. 10. (Color online) The top view of the same dynamical picture as in Fig. [9] 



destroys the soliton (its power at first increases, as it moves across the first PN barrier, and then decays to zero). 
At still higher fco, the soliton survives the kick, but in this case its trajectory is curvilinear in the plane of {X, Y)as 
shown in Fig. 1151 

It is relevant to compare the number of solitons in the patterns generated by the simulations for 9 = and 9 = tt/8. 
For both cases, these numbers are presented, as functions of projection ko x = fco cos # of the kick vector onto the 
X-axis, in Table [TTT1 It is seen that the dependences of the soliton number on kg x are quite similar for 9 = and 
9 = tt/8, barring the case of the destruction of the soliton (0 in the table), which occurs at 9 = tt/8, but does not 
happen for 9 = 0. 

Further, the evolution of the soliton's velocity for different strengths of the kick is shown in Fig. [16J Recall that, 
for 9 = tt/8, there are two domains of values of fco i n which the kicked soliton moves, separated by interval (|20j) 
where the soliton is destroyed by the kick. It is observed that the solitons accelerate and decelerate below and above 
the nonexistence interval ([2TJ]) . respectively, but eventually the velocity approaches a constant value. Moreover, the 
picture suggests that, as a result of the long evolution, the velocity is pulled to either of the two discrete values, w 2 
or » 3. As said above, these velocities (in other words, the values of fco which can directly produce such velocities, 
see straight horizontal lines in Fig. I16[) correspond to the single soliton moving along the X axis, rather than under 
an angle to it. The conclusion that the system relaxes to discrete values of the velocity is natural, as the dissipativc 
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FIG. 11. (Color online) (a): The field amplitude, |u (X, Z)\, in the cross section Y = 0, and (b) the evolution of the total 
power, for the kick's strength ko — 1.693 at 6 = 0, which leads to the establishment of the six-soliton pattern. 




(a) (b) 



FIG. 12. (Color online) (a) The distribution of the field amplitude, |w(-X", Z)\ in cross section Y = 0, and (b) the evolution of 
the velocity of the fundamental soliton kicked by ko = 2.1 at 6 — 0. 



system should give rise to a single or several isolated attractors, rather than a continuous family of states with an 
arbitrary velocity. 

A similar behavior is observed for other orientations of the kick, 9 = 7r/4 and 9 = 0, see the curves labeled by 
these values of 9 in Fig. [16j an d also Fig. 



12(b) 



Again, the velocity asymptotically approaches the same discrete 
values, close to 2 and 3, with the acceleration or deceleration below and above these values, respectively. In the case 
of 9 = 7r/4, in the established regime the free soliton moves in the diagonal direction. 



IV. COLLISIONS BETWEEN MOVING SOLITONS AND STANDING PATTERNS 



One of generic dynamical patters identified above features a standing multi-soliton structure and a freely moving 
soliton, which, due to the periodic boundary conditions, hits the standing structure from the opposite side, see Figs. 
[6j [9l [TTJ and [14] Two distinct scenarios of the ensuing interaction have been identified in this case, namely, the 
elastic collision, with the incident soliton effectively passing the quiescent structure (via a mechanism resembling the 
Newton's cradle, cf. Ref. (28j ) and reappearing with the original direction and velocity of the motion, and absorption 
of the free soliton by the structure, see Figs. 17(a) and |1 7(b)] respectively. 

More complex interaction scenarios were observed too, with several elastic or quasi-clastic collisions that end up 
with the eventual absorption of the free soliton. Outcomes of the collisions arc summarized in Table IIVI (the range 
of ko > 2.081 is not shown in the table, as only the single soliton exists in that case). At 1.766 < ko < 1.86, several 
clastic collisions, from 1 to 5, the number of which alternates in an apparently random fashion, are observed before 
the absorption is registered. This scenario is labeled " Newton's cradle with damping" in Table HVl At ko > 1.86, the 
collision is elastic and persists to occur periodically, as in the case of the ordinary Newton's cradle (so named too in 
Table [TV]), see Fig. QB 

A special case is the one corresponding to the creation of the largest number of solitons at ko = 1.693, as shown 



above. The collision pattern is quite complex in this case, as shown in Fig. 19(a) Both clastic collisions (at Z ~ 50 
and 100) and absorptions (at Z ~ 238) arc observed. An unexpected feature of the process is the reversal of the 
direction of motion of two solitons around Z ~ 150. The whole patterns eventually relaxes into an array built of six 
solitons, which is confirmed by the power-evolution plot in Fig. |19(b)[ 

In the case of elastic collisions between two solitons recurring indefinitely long, we have checked if the velocity of 
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FIG. 13. (Color online) The distribution of the field amplitude, \u (X, Y)\ (a) and phase (b) in the pattern produced by kicking 
the fundamental soliton in the diagonal direction (9 — 7r/4) very hard, with ko = 100. The picture corresponds to evolution 
distance Z = 49.98. 




FIG. 14. (Color online) The local field amplitude, \u(X,Z)\, in the cross section Y = 0, for ko = 1.9743 and 6 — tt/8. 



the transmitted soliton is the same as that of the incident one. To this end, it is relevant to consider the case of 
fco = 1.974 and 6 = 0. After the first round trip of the emitted soliton, the one original is quiescent while the emitted 
one is running into it. It is seen in Fig. 20(a) that each collision leads to the exchange of velocities between the two 
solitons, as for colliding hard particles. In addition, we identified the velocity of the center of mass of the two-solitons 
set. Figure [20 (b ) | shows that the latter velocity gradually increases within the sequence of collisions, approaching one 
of the above-mentioned discrete values characteristic to the established motion of single solitons. In the present case, 
the asymptotic velocity is close to 1. 



V. CONCLUSIONS 



The subject of this work is the mobility of 2D dissipativc solitons in the CGL (complex Ginzburg-Landau) equation 
which includes the spatially periodic potential. This equation models bulk lasing media with built-in transverse 
gratings. The soliton was set in motion by the application of the kick, which corresponds to a tilt of the seed beam [29| . 
The mobility implies the possibility to generate oblique laser beams in the medium. Further, the advancement of the 
kicked soliton may be used for the controllable creation of various arrayed patterns in the wake of the soliton hopping 
between the potential cells. The depinning threshold, i.e., the smallest strength of the kick which sets the quiescent 
soliton into motion, was found by means of simulations, and also with the help of the analytical approximation, 
based on the estimate of the condition for the passage of the kicked object across the PN (Peierls-Nabarro) potential. 
The dependence of the threshold on the orientation of the kick with respect to the underlying lattice potential was 
studied too. Various pattern-formation scenarios have been identified above the threshold, with the number of solitons 
in stationary arrayed patterns varying from one to six. Freely moving solitons may eventually assume two distinct 
values of the velocity, which represent coexisting attractors in this dissipative system. Also studied were elastic and 
inelastic collisions between the free soliton and stationary multi-soliton structures, with two generic outcomes: the 
quasi-clastic passage, like in the case of the Newton's cradle, and absorption of the free soliton by the quiescent 
structure (sometimes, after several passages). A natural extension of the present work may deal with the dynamics 
initiated by the application of the kick to vortices pinned by the underlying grating. 
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FIG. 15. (Color online) Examples of curvilinear trajectories of the soliton for 9 = n/S and two large values of the kick, fco = 6.05 
and 10. For the references, straight lines with sloped 9 — tt/8 and 9 = 7r/4 are displayed too. 



Range of fco 


Range of the 


Number of solitons 


Number of solitons 




projection ko x 


= 


6 = tt/8 


[0, 1.816] 


[0, 1.6778] 


1 


1 


1.974 


1.8237 


3 


3 


2.1 


1.9401 


2 


2 


[2.224, 4.569] 


[2.0547,4.2212] 


1 


1 


[4.816, 5.804] 


[4.4494, 5.3622] 


1 






TABLE III. The number of solitons versus the ko x component of the kick vector, for orientations 8 = and tt/8. 
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FIG. 16. (Color online) The soliton's velocity as a function of Z at various values of fcrj and 6. 
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TABLE IV. Collision types versus the magnitude of the initial kick, fco. 
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FIG. 18. (Color online) An example of periodic elastic collisions according to the Newton's-cradle scenario, at fco = 1.867 and 
9 = 0. 
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FIG. 19. (Color online) (a): An example of the complex interaction of the free soliton with the standing structure at ko = 1.693 
and 9 — 0. (a) The distribution of the field, \u (X, Z)\, in the cross section Y = 0; (b) the total power versus Z. In panel (b), 
horizontal lines show the multiple powers of the quiescent fundamental soliton. 




FIG. 20. (Color online) Velocities of each soliton (a) and center of mass (b) versus Z for the pair of periodically colliding 
fundamental solitons. Here ko = 1.974 and 9 — 0. 



